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Abstract 

Interactions within the cosmic medium modify its equation of state. We discuss implications 
of interacting dark energy models both for the spatially homogenous background and for the 
perturbation dynamics. 
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I. INTRODUCTION 



According to our current understanding the presently observable Universe is dynamically 
dominated by two so far unknown substances, Dark Matter (DM) and Dark Energy (DE). Of 
major interest is a precise knowledge of the DE equation of state (EOS) including its possible 
time dependence. A definitely time dependent EOS would rule out a "true" cosmological 
constant which still is the most favored DE candidate. A related problem is to understand 
why the energy densities of both components are of the same order of magnitude today (the 
"coincidence problem" ) . 

Most approaches in the field rely on an independent evolution of DE and DM. Given the 
unknown nature of both DE and DM one may argue, however, that an entirely independent 
behavior is a very special case At least, there is no compelling reason to exclude 

interactions from the outset. Unified models such as the Chaplygin gas models even try to 
understand DE and DM as different manifestations of one single fluid jjj. 

Our aim here is to point out the role of interactions within the dark sector for the 
cosmological dynamics. As far as the coincidence problem is concerned, an interaction may 
naturally lead to a fixed ratio of the energy densities of DM and DE. An interaction will also 
influence the perturbation dynamics. In particular, it modifies the adiabatic sound speed 
and the large scale perturbation behavior. It affects the time dependence of the gravitational 
potential and hence the Integrated Sachs- Wolfe effect (ISW) which is observable through 
the lowest multipoles of the anisotropy spectrum of the cosmic microwave background. Fur- 
thermore, a coupling between DM and DE is expected to be relevant for holographically 
determined DE and for the big rip scenario. 

II. INTERACTING DARK ENERGY 

An interaction between DM and DE is most conveniently modelled as a decay of the 
latter into the former, 

p M + 3Hp M = Tp x , px + 3H (1 + w x ) px = ~^px , (1) 

where pu and px are the energy densities of DM and DE, respectively, H is the Hubble 
expansion rate and wx is the EOS parameter of the DE. The (not necessarily constant) 
decay rate is denoted by Y. Einstein's field equations for spatially flat FRLW cosmologies 
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where p = pu + Px is the total energy density. The quantity H is related to the deceleration 
parameter q by q = — 1 — -jp. It is obvious, that neither the Hubble parameter H nor the 
deceleration parameter is directly affected by the interaction (there is an indirect influence, 
however, since the coupling modifies the ratio — ). It is only the second derivative of the 

Px 

Hubble parameter, equivalent to the third derivative of the scale factor a (where H = -) in 
which the decay rate enters explicitly, 
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This corresponds to the fact that T enters the luminosity distance dt, = (1 + z) J ^ only 
in third order of the redshift z (cf. Refs. ^,0]). Consequently, the interaction becomes 
observable if the cosmological dynamics is not only described by the present values H and 
go of H and q, respectively, but additionally by parameters that involve at least Hq, the 
present value of H. The recently introduced "statefinder" parameters (f| [?J are of this type 
and represent useful tools to characterize interacting models 



III. CONSTANT DENSITY RATIO 

Assuming wx = const there exists a special solution of the system (JIJ for which k = 
= kq = const. This solution which is of interest for the coincidence problem |2] requires 
a decay rate 

r = -3H-^-w x . (4) 

The EOS parameter wx has to be negative to obtain T > 0, i.e., a decay of the component 
X into the matter component. Under this condition the energy densities scale as 

p M , p x , p oc a-^ 1+w \ w = . (5) 

Here, w is the overall EOS which, compared with wx, is reduced in magnitude by a factor 

2 

1 + k . The time dependence of the scale factor in the spatially flat case is a oc i 3 < 1+ ™) . 

In a next step we realize that the chain of arguments leading to (0) (via (0) and (@J)) may 
be reversed. The resulting statement is that any cosmic EOS w can be interpreted as an 



effective EOS of an interacting mixture of a component with EOS wx and a dust component 
(a generalization to an arbitrary constant EOS is possible). With other words, for any 
one-component description of the cosmic medium there exists an equivalent interacting two- 
component model which generates the same background dynamics. This equivalence (which 
implies a degeneracy with respect to observational results) will turn out to be useful on the 
perturbative level. It provides us with an internal structure for any one-component cosmic 
medium and represents a systematic way to generate non-adiabatic perturbations which we 
shall consider later on. 

An example of how a solution kq = const may be approached as the result of a dynamical 
evolution of the energy density ratio k can be sketched as follows [jj. For a decay rate 
T = 3Hb 2 (1 + k) with a parameter b 2 = const, the set of equations (Q) admits two stationary 
solutions, k + and k_, for the ratio k of the energy densities. These have the properties 
> and k + k_ = 1. For k = k + = k_ we have pu > Px, i-e., matter dominance, 
while for k, = k~ the reverse relation Pm < Px holds, equivalent to a DE dominated phase. 
There exists a dynamical solution k it) according to which the density ratio evolves from 
k = kZ 1 for small values of the scale factor (matter dominance) towards a stable, stationary 
solution k = k_ for large values of the scale factor (DE dominance). The deceleration 
parameter changes from positive values at hT} to negative values at k_. This means, the 
interaction drives the transition from a matter dominated era to a subsequent phase with 
accelerated expansion. The solution k = k_ can be identified with the solution kq, given by 
(jHJ) for b 2 = —kqWx (1 + Ko) 2 - Consequently, a stationary ratio for a density ratio which 
generates accelerated expansion is obtained as the result of a dynamical evolution. 

As already mentioned, the overall equation of state w in (jSJ) is reduced in magnitude 
compared with wx- In particular, one may ask whether it is possible to have w > — 1 while 
wx < — 1. To this purpose it is useful to consider the basic equations (JTJ) for the simple case 
T = const. Even if at the start of the decay the evolution is dominated by px, after a short 
time pm will be approximately 
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The calculations leading to this result are similar to those at the end of an out-of-equilibrium 
decay of (e ff e Cti ve ly ) d U8t mat te r into Nation ta inflation — Q Q. It is 



obvious, that the ratio ^ in © can only be constant for H = const, i.e., for exponential 



expansion. Combining now (JUJ) with (j^J) we obtain a direct relation between wx and /to, 

wx = — • (7) 
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This configuration has wx < — 1 necessarily. It approaches —1 only in the limit of a vanishing 
matter component. While a non-interacting DE component with Wx < — 1 is known to lead 



to a singularity in a finite time 
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13j (see, however, Ref. jlj|), our example demonstrates 
that a suitable interaction makes the expansion de Sitter like, thus avoiding the big rip. 

m 



Similar conclusions for interacting phantom energy have been obtained in Refs. 

(Alternative settings with interactions that lead to a transition from wx > — 1 to Wx < — 1 

f| f] 

were discussed in Refs. |18j and |19|). 

It is also interesting to see how a solution (JHJ) with k = /to is related to the holographic 
bound on the DE Q. With an infrared cutoff scale L the DE density is p 1 ^ 1 = 3c 2 MpL~ 2 
where Mp = 8itG (cf. Ref. jsi]]). It seems natural to identify L with the Hubble scale H -1 . 
With this choice it follows from Friedmann's equation that both px and pu are proportional 
to H 2 . While this property was recently used to discard a cutoff L = H~ l for independently 



evolving components 



models 
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2 it ] . the situation is entirely different in the context of interacting 



23| . Namely, the same dependence of both px and p m on H can be regarded as 



the result of a suitable interaction between these components in the sense described by (@J) 



and (0) • With a cutoff L = H 1 we have k 



pm 
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k = const from the outset . Moreover, it 



immediately relates the parameter c in the expression for p£ to the ratio Kn by n 
It follows that c 2 < 1, different from the case of non-interacting models |2l| where L is 
identified with the future event horizon (for an interacting model in the latter context see 
Ref. The condition for accelerated expansion, — 1 < w < —1/3, then translates into 

- (/t + 1) < w x < ~ (/«o + 1) equivalent to - < w x < ■ (8) 

We emphasize that the interaction is essential here to have acceleration. There is no non- 
interacting limit. Without interaction acceleration is impossible, i.e., acceleration is a pure 
interaction phenomenon in this picture. 



IV. PERTURBATIONS 

Adiabatic perturbations of the cosmic fluid are characterized by p = jp, where p and p 
denote the total pressure and energy density perturbations, respectively. Both quantities 
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are related by the adiabatic sound speed If the substratum is understood as composed 
of two subcomponents with component X decaying into component M, the decay rate T 
may fluctuate. Fluctuating decay rates are known to produce curvature perturbations in 



certain inflationary scenarios [25|, |2y, l27| . Here, in a different context, they are used to 
characterize internal structure in the dark sector of the cosmic substratum For the 

medium discussed here we find 

where (jjj) denotes fluctuations about the (constant) background ratio A fluctuating 
fractional decay rate ^ necessarily modifies the adiabatic sound speed. In particular, for a 
fluctuation of the type (^) = A^ we obtain an effective sound speed square Cg ff = j + A. 
Even if ? is negative, a non-adiabatic contribution due to an internal structure may give rise 
to a physically acceptable effective sound speed. On the perturbative level the cosmic EOS 
may be different from p = wp. Fluctuations of the ratio ^ can be regarded as fluctuations of 
the EOS which give rise to pressure perturbations p = wp + wp, where w = (^) • Thus, the 
equivalent two-component picture represents a simple method to construct variations of the 
cosmic EOS while keeping % = and wx = const. On this basis one may reconsider fluid 
models of the cosmic substratum which occasionally are abandoned since a purely adiabatic 
treatment does not correctly reflect the observational situation. Taking into account an 
internal structure of the medium, which is equivalent to go beyond the limits of an adiabatic 
analysis, will generally result in a different picture. 

The internal structure, modelled here through a fluctuating decay rate, will also influence 
the large scale perturbation behavior. Large scale perturbations are conveniently described 
by the quantity C which represents a curvature perturbation on hypersurfaces of constant 
. On large perturbation scales ( obeys the equation (cf. Refs. j^] - [3^ ) 



energy density 



p — £p 

C = -n— p - ■ (10) 
p + p 

For adiabatic perturbations p = jp the quantity ( is conserved. As soon as there are non- 
adiabatic perturbations, ( will no longer be constant. For our model with w > — 1 (and, for 
simplicity, time independent perturbations w = (j^) ) we find 

f = c,-^(4W-V (id 



1 + w \ 3H I \a 



A fluctuating fractional decay rate gives rise to a change in the curvature perturbation 
which is logarithmic in the scale factor. The subscript i denotes some initial value. The 
circumstance that ( is time dependent has consequences for the time dependence of the 
(gauge-invariant) gravitational potential for which we obtain 



5 + 3w " 25 + 3w \ 3H 



In I- 
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While the first constant term on the right-hand side is a standard result for adiabatic per- 
turbations [2^, the fluctuating decay rate introduces a logarithmic dependence on the scale 
factor. The time dependence of $, 

* = £_!_ Cn'a, (13) 

25 + 3wV©y 

determines the ISW. Different from the pure adiabatic case where $ = there is an ISW 
in this case. It is interesting to compare the result (|T2|) with the corresponding expression 
of the ACDM model. In the long time limit one finds for the latter $acdm ~ ~\d^f^^ 
where is some initial value and pa = const. This means, <3?acdm shows a stronger 
time dependence than $ in f!12|) . Apparently, our interacting DE model predicts a smaller 
ISW compared with the non-interacting ACDM model. A smaller ISW is equivalent to a 
suppression (compared with the prediction of the ACDM model) of the lowest multipoles 
in the anisotropy spectrum of the cosmic microwave background. Such kind of suppression 
is required by the data and different mechanisms have been invented to account for this 
phenomenon 

Q Q, H Q- 

Here we argue that it may well be the consequence of an 

interaction in the dark sector. 



V. SUMMARY 

Interacting DE models are more general than non-interacting ones. A coupling between 
DE and DM enriches the cosmological dynamics. It reveals potential degeneracies in inter- 
preting the observational data and provides a method to refine the equation of state of the 
cosmic substratum. An interaction is particularly useful to address the coincidence problem. 
It also opens new perspectives on holographic DE and on the big rip scenario. Furthermore, 
it introduces a non-adiabatic feature into the perturbation dynamics with consequences for 
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the effective sound speed and the time dependence of the gravitational potential. The in- 
teraction should be detectable by suitably refined luminosity distance measurements and by 
certain features in the anisotropy spectrum of the cosmic microwave background. 
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